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Abstract. We report on the detection of three Satum-mass planets discovered with the HARPS instrument. HD 93083 shows 
radial-velocity (RV) variations best explained by the presence of a companion of 0.37 Afj„p orbiting in 143.6 days. HD 101930 b 
has an orbital period of 70.5 days and a minimum mass of 0.30 Mj^p. For HD 102117, we present the indep endent detect i on of a 
companion with m 2 sini = 0.14 Mj„p and orbital period P = 20.7 days. This planet was recently detected bv lTinnev et al.U2004) . 
Activity and bisector indicators exclude any significant RV perturbations of stellar origin, reinforcing the planetary interpreta¬ 
tion of the RV variations. The radial-velocity residuals around the Keplerian fits are 2.0, 1.8 and 0.9 m s“* respectively, showing 
the unprecedented RV accuracy achieved with HARPS. A sample of stable stars observed with HARPS is also presented to 
illustrate the long-term precision of the instrument. All three stars are metal-rich, confirming the now well-established relation 
between planet occurrence and metallicity. The new planets are all in the Satum-mass range, orbiting at moderate distance 
from their parent star, thereby occupying an area of the parameter space which seems difficult to populate according to planet 
formation theories. A systematic exploration of these regions will provide new constraints on formation scenarios in the near 
future. 

Key words, stars: individual: HD 93083, HD 101930, HD 102117 - stars: planetary systems - techniques: radial velocities - 
techniques: spectroscopic 


1. Introduction 

The HARPS instrument is the new ESO high-resolution (R 
- 115,000) fiber-fed echelle spectrograph, mainly dedicated 
to planet search programmes and asteroseismology. It has al¬ 
ready proved to be the most precise spectro-velocimeter to 
date, reaching an instrumental RV accuracy of ~ 1 m s“' 

llMavor et '^l2003t Ip^e et al.ll2004t Is antos et aDl2004al) . and 

even better on a short-term basis. This opens a new field in the 
search for extrasolar planets, allowing the detection of com¬ 
panions of a few Earth masses around solar-type stars. Indeed, 
the increase of the planet frequency towards very low masses is 
confirmed by the recent discovery of Neptune-mass exoplanet s 
llMcArthur et al.ll2004l ISantos et al.ll20043lButler et al.l200^ . 


Send offprint requests to: C. Lovis 
* Based on observations made with the HARPS instrument on the 
ESO 3.6 m telescope at La Silla Observatory under programme ID 
072.C-0488(E). 


Moreover, the combination of CORALIE dOueloz et alJl200ol) 
and HARPS data, extending over several years, will allow the 
characterization of long-period planets and complex planetary 
systems. The low-mass and long-period regions of the exo¬ 
planet parameter space will therefore be under close scrutiny in 
the coming months and years. This will improve our knowledge 
of the planet distribution in the mass-period diagram and will 
all ow comparisons with theoretical predictions (see for exam ¬ 

ple lArmitage et al.l2002llida & Linl2004rAlibert et a 1 1200.51) . 

In this paper we present the discovery of three plane¬ 
tary companions to the stars HD 93083, HD 101930 and 
HD 102117. Interesting characteristics of these planets include 
a mass in the Saturn-mass regime and below, and an orbital dis¬ 
tance to the star o f 0.1-0.5 AU. Ac cording to recent planet for¬ 
mation scenarios jlda & Linl2004 . they are situated in a region 
of the mass-period diagram that seems difficult to populate. 
The systematic exploration of this mass and distance regime is 
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therefore an important test for the core-accretion and migration 
theories of planet formation. 

This paper is structured as follows. Section 2 describes 
the observations, the data reduction process and discusses the 
long-term precision of the HARPS measurements. The physi¬ 
cal properties of the parent stars are presented in Sect. 3, fol¬ 
lowed by the RV measurements and orbital solutions in Sect. 4. 
We discuss in the last section the characteristics of these new 
planets in the context of the already-known properties of exo¬ 
planets and highlight some theoretical questions that might be 
answered by radial-velocity surveys in the near future. 

2. Observations and measurement precision 

HD 93083, HD 101930 and HD 102117 are all members of 
the HARPS high-precision RV sample. The aim of this sur¬ 
vey is to obtain RV measurements with photon errors below 
1 m s“' to detect very low-mass extrasolar planets. The stars 
in this sample have been selected from the CORALIE planet 
search database for being non-evolved, having low projected 
rotational velocity (v sin i < 3 km s“') and exhibiting low activ¬ 
ity levels (log < -4.7). These criteria should eliminate most 
stars showing large intrinsic RV variations and select only very 
quiet, chromospherically inactive stars, so that stellar RV jit¬ 
ter does not hide possible planetary signals down to 1-2 m s '. 
Obviously, the photon noise error on the radial velocity must 
also remain below 1 m s '. For each HARPS spectrum we 
compute the photo n-limited , ultirn ate R V precision u s ing th e 
formulae given by IConnesI (Il985l) and iBouchv et alJ ll200ll) . 
HARPS typically delivers 0.5-1 m s ' photon-limited accuracy 
at a S/N ratio of ~100 per pixel (0.8 km s“') at 550 nm on solar- 
type stars. The exact number mainly depends on the depth of 
absorption lines and the projected rotational velocity. 

The radial velocities for HD 93083, HD 101930 and 
HD 102117 have been obtained with the standard HARPS re¬ 
duction pipeline, based on the cross-correlation with a stellar 
template, the precise nightly wavelength calibration with ThAr 
spectra and the tracking of instrumental drift s with the simulta¬ 
neous ThAr technique (iBaranne et al.lll9^ . Of particular in¬ 
terest is the fact that the nightly instrumental drifts always re¬ 
main below 1ms' for HARPS due to the high instrumental 
stability. 

To assess the long-term precision of the instrument, we 
selected a sample of 12 stable stars, which have been fol¬ 
lowed regularly over the past two years. These stars include 
HD 55 (K5V, y=8.49), HD 1581 (F9V, y=4.23), HD 6673 
(KIV, y=8.84), HD 20794 (G8V, y=4.26), HD 28471 
(G5V, y=7.89), HD 30278 (G8V, y=7.61), HD 44594 (G4V, 
y=6.61), HD 82342 (K3/K4V, y=8.31), HD 114853 (G2V, 
y=6.93), HD 162396 (F8V, y=6.19), HD 196761 (G8/K0V, 
y=6.36) and HD 215152 (KOV, y=8.11). For each of these 
stars, several measurements have been gathered, spanning at 
least 300 days, and the mean RV was subtracted to enable a 
comparison between different stars. The first data points were 
acquired during the commissioning periods of the instrument 
in February and June 2003. Between then and the more recent 
measurements, many modifications have been made to the in¬ 
strument, especially the exchange of the ThAr reference lamps. 



52600 52800 53000 53200 53400 

JD - 2400000 


Fig. 1. Radial velocities for a sample of 12 stable stars from 
the HARPS high-precision programme, observed over about 
700 days. For each star, the mean RV has been subtracted. 
The overall RMS of 1.81 m s ' is made of different contri¬ 
butions, among which photon noise (~0.8 m s“'), calibration 
noise (~0.8 m s“'), guiding errors (0.2-0.6 m s“') and stellar 
acoustic modes (0.5-2.0 m s“'). Seasonal averages have been 
computed (bottom part). The absence of long-term drifts and 
the dispersion of only 0.64 m s“' illustrate the long-term sta¬ 
bility of ThAr reference lamps. 

In spite of that, the global RV dispersion of the measurements 
amounts to 1.81 m s“' (see Fig0. This includes many different 
effects, which can be estimated as follows. 

- The mean photo n error, computed according to 
iBouchv et al.lll200il) . is 0.8 m s '. 

- The zero-point wavelength calibration, given by nightly 
ThAr reference spectra, is presently determined with an ac¬ 
curacy of ~0.8 m s“'. 

- Guiding errors are responsible for 0.2-0.6 m s ' of extra 
dispersion, due to the non-optimal tuning of the guiding 
software. This issue has now been solved and the guiding 
noise should remain below 0.3 m s ' in the future. 

- The largest contribution to the global RV dispersion comes 
from the star itself, mainly due to acoustic modes (p- 
modes). For our 12 stars we estimate this contribution to 
be 0.5-2.0 m s ', depending on spectral type and evolu - 
tionary status tsee lMavor et al.ll2003URouchv et al]l200.5[) . 
The observational strategy has been recently optimized to 
minimize these stellar oscillations, essentially by integrat¬ 
ing over long enough periods to cover more than 1 -2 oscil¬ 
lations (3-15 minutes, depending on spectral type). 

Adding all these effects quadratically, we obtain a total dis¬ 
persion that is very close to the measured value of 1.81 m s“' 
(the exact number depending on the adopted mean value for 
the oscillation noise). This leaves only little space (probably 
less than 1 m s“') for other star-related variations, such as the 
presence of undetected planetary companions, star spots and 
activity-induced jitter. 
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Table 1. Observed and inferred stellar parameters for the 
planet-hosting stars presented in this paper. 


Parameter 

HD 93083 

HD 101930 

HD 102117 

Sp 


K3V 

KIV 

G6V 

V 

[mag] 

8.30 

8.21 

7.47 

B-V 

[mag] 

0.945 

0.908 

0.721 

n 

[mas] 

34.60 

32.79 

23.81 

My 

[mag] 

6.00 

5.79 

4.35 

Teff 

[K] 

4995 ± 50 

5079 ± 62 

5672 ± 22 

logg 

[cgs] 

4.26 ± 0.19 

4.24 ±0.16 

4.27 ± 0.07 

[Fe/H] 

[dex] 

0.15 ±0.06 

0.17 ±0.06 

0.30 ± 0.03 

L 

[Lo] 

0.41 

0.49 

1.57 

M. 

[Mo] 

0.70 ± 0.04 

0.74 ± 0.05 

1.03 ±0.05 

V sin i 

[km s“'] 

0.9 

0.7 

1.5 

logRjlK 


-5.02 ± 0.02 

-4.99 ± 0.02 

-5.03 ± 0.02 

-^rot 

[days] 

48 

46 

34 


Among the above-mentioned error sources, three of them 
should be of less importance in the future, permitting an even 
better precision. Firstly, guiding errors are no longer a limiting 
factor. Secondly, the new observing strategy should allow us 
to maintain the stellar oscillation ’’noise” below 0.5 m s ' for 
K-dwarfs and 1.5 m s * for slightly evolved, early G-dwarfs 
(worst case). Indeed, comparisons between measurements on 
a given star taken before and after the strategy change have 
already shown a signihcant reduction of the RV dispersion. 
Finally, improvements in the wavelength calibration are still 
ongoing, in particular the elaboration of a more accurate list of 
ThAr laboratory wavelengths, which will lead to a lower cali¬ 
bration noise. The improved wavelength solutions will then be 
applied to all data (past and future), ensuring an even better 
night-to-night stability. 

To study the behavior of reference ThAr spectra over time, 
we have computed seasonal averages (see FiglU bottom). The 
dispersion of these seasonal mean values amounts to 0.64 
m s ', showing that ThAr lamps are able to deliver an absolute 
wavelength reference at that level of precision over long peri¬ 
ods of time, even in the case of lamp exchange or lamp current 
variations (both occurred during the hrst year of HARPS op¬ 
eration). In conclusion, we are therefore very conhdent of the 
HARPS ability to maintain a long-term accuracy at the level of 
1 m s“'. 

3. Parent star characteristics 

The basic properties of HD 93083 (K3V, y=8.30), HD 101930 
(KIV, y=8.21)andHm02117(G6V, y=7.47) are given by the 
Hipparcos catalogue ilFsall997h . Physical properties have bee n 
derived following the method by ISantos et~ ] il200lH2004l3) . 
using a standard local thermodynamical equilibrium (LTE) 
analysis. This study provides precise values for the effective 
temperatures, metallicities and surface gravity. 

From the colour index, the measured effective temperature 
and the corresponding bolometric correction, we estimate the 
star luminosities and we then interpolate the masses and ages in 
the grid of Geneva stellar evolutionary rnodels with appropriat e 
metal abundances dSchaller et alJll992t ISchaerer et al.lll 993l) . 



Wavelength [A] 

Fig. 2. Call H (T = 3968.47 A) absorption line region for 
HD 93083 (top), HD 101930 (middle) and HD 102117 (bot¬ 
tom). Reemission peaks are absent or very weak in these three 
stars, indicating low chromospheric activity. 

The projected rotational velocity v sin i is also computed us¬ 
ing the c alibration of t he CO RALIE cross-correlation function 
given bv ISantos et al.l(l2002h . Tabled gathers the photometric, 
astrometric, spectroscopic information and inferred quantities 
for HD 93083, HD 101930 and HD 102117. 

The results show that the masses of these stars range from 
0.7 to 1.03 Mq. Interestingly, they all have a metallicity higher 
than solar ([Ee/H] = 0.15, 0.17 and 0.30 respectively), con- 
hrming the now well-es tablished relation be tween planet oc¬ 
currence and metallicity dSantos et al.l2004bl) . 

We also compute from the spectra an activity indicator 
Sharps measuring the re-emission flux in the Call H and K 
lines. This indicator has be en calibrated on the Mount Wilson 
scale jBaliunas et alJll993) using stars with well-known activ¬ 
ity leve ls. We then correct f or the photospheric contribution fol¬ 
lowing ^o^e^^O ill984li and obtain the well-known activity 
index logRjjj.. This index represents a useful tool to estimate 
the stellar RV jitter expected for each star due to rotational 
modulat ion of star spots or oth er active regions on the stellar 
surface jSaar & Donahueil997h . Erom this indicator we derive 
an estimati on of the stellar rot ation period following the cali¬ 
brations bv iNoves et al.l ilT984l) . Although these relations have 
been shown not to be reliable for ages greater than ~2 Gyr 
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jPace & Pasauinl2004l) . they are still helpful in distinguishing 
between active, young stars and chromospherically quiet, old 
stars. 

All three stars examined in this paper exhibit low chro¬ 
mospheric activity (log/?[jj^ = -5.0, see FigEJ. This is further 
supported by the non-detection of X-ray emission from these 
stars, despite their relatively small distance. Together with the 
measured small vsin/ (< 1.5 km s“^), this indicates very low 
activity-induced RV jitter (< 1-2 m s“*). Moreover, the analy- 
sis of the bisector s hape of the cross-correlation function (see 
lOueloz et al.l I200l1l shows no variations in the CCF profile 
down to the photon noise level, giving strong support to the 
planetary interpretation of the RV variations. 


4. Radial-velocity data and orbital solutions 

4.1. HD 93083 

Sixteen radial-velocity measurements have been obtained for 
HD 93083 spanning about one year. Exposure time was 15 
min on average, yielding a typical S/N ratio of 130 per pixel at 
550 nm. The mean photon noise error on a single data point is 
0.4 m s"', to which we quadratically added the corresponding 
calibration error (see Sect. 2). The list of all RV measurements 
for HD 93083 can be found in Tabled Fig|3shows the radial 
velocities folded to a period P = 143.58 days, given by the best 
Keplerian fit to the data. The RV semi-amplitude is k = 18.3 
m s“' and the eccentricity e - 0.14. These parameters lead to a 
minimum mass m 2 sin i = 0.37 Mjup for the planet. All relevant 
data for this planet are listed in TableEl 

The weighted rms of the residuals around the fit is 2.0 
ms*, which is larger than the internal errors. More data points 
are needed to establish whether this extra dispersion is intrinsic 
to the star or whether it could be explained by the presence of 
a third body in the system. The latter hypothesis is reinforced 
by the slow RV drift that seems to be present in the residuals 
of the radial velocities around the orbit, as can be seen in Fig|3 
(bottom). 

4.2. HD 101930 

We have obtained 16 HARPS radial-velocity measurements of 
HD 101930 over a time span of ~ 1 year. The exposure time was 
15 min, yielding a typical S/N ratio of 110 per pixel at 550 nm. 
The mean photon noise error on a single data point is 0.5 ms"*, 
to which we quadratically added the calibration error (see Sect. 
2). Tablel^gives the list of all RV measurements obtained for 
HD 101930. 

The RV variations are best explained by a Keplerian orbit 
with period P - 70.46 days, semi-amplitude k - 18.1 m s"* 
and eccentricity e = 0.11. From these orbital parameters and 
the primary mass given in Table [2 (0.74 Mq), we derive a 
minimum mass m 2 sin / = 0.30 Mjup for the companion. The 
residuals around the fit have a weighted rms of 1.8 m s"*, 
slightly larger than internal errors. The extra dispersion might 
be caused by stellar RV jitter. Future measurements will al¬ 
low us to investigate this point more precisely. Fig@]shows the 
phased radial velocities for HD 101930 together with the resid- 


Table 3. Radial-velocity values and error bars for HD 93083. 
All data are relative to the solar system barycenter. 


JD-2400000 

RV 

[km s"*] 

Uncertainty 
[km s"*] 

53017.84496 

43.6470 

0.0014 

53036.78231 

43.6607 

0.0013 

53046.66898 

43.6638 

0.0013 

53049.71044 

43.6651 

0.0013 

53055.74405 

43.6588 

0.0013 

53056.75374 

43.6587 

0.0014 

53060.75307 

43.6586 

0.0013 

53145.55149 

43.6324 

0.0010 

53151.57396 

43.6345 

0.0010 

53202.47605 

43.6581 

0.0009 

53205.45320 

43.6594 

0.0009 

53310.85564 

43.6478 

0.0009 

53343.85438 

43.6551 

0.0008 

53375.83553 

43.6369 

0.0008 

53377.84184 

43.6362 

0.0008 

53400.79648 

43.6266 

0.0012 


HD 93083 HARPS 



7 0.004 

m 

g 0.002 

A 0 

U -0.002 

^ -0.004 

53000 53100 53200 53300 53400 

JD-2400000 [days] 



Fig. 3. Phased radial velocities for HD 93083. The best 
Keplerian fit to the data gives a minimum mass of 0.37 Mjup 
and an orbital period of 143.6 days for the planet. 


uals around the orbit, and TableElcontains all relevant parame¬ 
ters for HD 101930 b. 

4.3. HD 102117 

For this star we gathered 13 radial-velocity measurements with 
HARPS spanning about 380 days between January 2004 and 
January 2005. These are listed in Table |3 The spectra, ac¬ 
quired in 10-15 min exposures, have a mean S/N ratio of 120 
per pixel, corresponding to a photon error of 0.5 ms* on this 
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Table 2. Orbital and physical parameters for the planets presented in this paper. 


Parameter 

HD 93083 b 

HD 101930 b 

HD 102117 b 

P 

[days] 

143.58 ± 0.60 

70.46 ±0.18 

20.67 ± 0.04 

T 

[JD-2400000] 

53181.7 ±3.0 

53145.0 ±2.0 

53100.1 ±0.1 

e 


0.14 ±0.03 

0.11 ± 0.02 

0.00 (±0.07) 

V 

[km s“'] 

43.6418 ± 0.0004 

18.3629 ± 0.0003 

49.5834 ± 0.0003 

(Jj 

[deg] 

333.5 ± 7.9 

251 ± 11 

162.8 ± 3.0 

K 

[m s“‘] 

18.3 ±0.5 

18.1 ±0.4 

10.2 ± 0.4 

a\ sin i 

1 

O 

0.239 

0.116 

0.019 

f(m) 

[lO-^Me] 

0.088 

0.042 

0.0023 

m 2 sin i 

[Afjup] 

0.37 

0.30 

0.14 

a 

[AU] 

0.477 

0.302 

0.149 

•^meas 


16 

16 

13 

S pan 

[days] 

383 

362 

383 

a (0-C) 

[m s“'] 

2.0 

1.8 

0.9 


Table 4. Radial-velocity values and error bars for HD 101930. 
All data are relative to the solar system barycenter. 


JD-2400000 

RV 

[km s '] 

Uncertainty 
[km s'] 

53038.78552 

18.3687 

0.0014 

53053.89075 

18.3485 

0.0014 

53055.83816 

18.3478 

0.0013 

53145.58216 

18.3582 

0.0010 

53147.55920 

18.3617 

0.0009 

53150.66830 

18.3660 

0.0009 

53154.58438 

18.3716 

0.0010 

53159.54707 

18.3811 

0.0009 

53201.47066 

18.3474 

0.0009 

53203.47312 

18.3449 

0.0009 

53206.46791 

18.3452 

0.0009 

53217.51591 

18.3596 

0.0009 

53218.51090 

18.3641 

0.0009 

53310.86794 

18.3776 

0.0009 

53344.86147 

18.3403 

0.0009 

53400.81929 

18.3566 

0.0010 


HD 10t930 


HARPS 




Star. Note that the photon noise at a given S/N ratio is the same 
as for HD 101930, although HD 102117 has an earlier spec¬ 
tral type (G6 vs. Kl), and therefore less numerous and broader 
spectral lines. This is probably due to the higher metallicity of 
HD 102117, leading to deeper absorption lines. After quadrat- 
ically adding the calibration error, we obtain a total error of 
~0.9 m s ' for each data point. 

The best-ht Keplerian orbit explaining the RV variations 
has a period P - 20.67 days, a semi-amplitude k - 10.2 m s“* 
and an eccentricity e equal to 0. We derive a companion min¬ 
imum mass m 2 sin / = 0.14 Mjup, which is one of the lightest 
extrasolar planets known to date. The weighted rms around 
the fit is only 0.9 m s“*, equal to internal errors. This shows 
that intrinsic stellar effects, such as activity-related jitter, are 
negligible for this star at a level of a few ten s of cm s~'. 
HD 10 2117 b has been previously announced bv lTinnev et al.l 
ll2004 . who find orbital parameters very similar to ours, al¬ 
though with slightly larger uncertainties due to the larger error 
bars ( (TQ-c = 3.3 m s“* ). The combination of both data sets 
show that the radial velocities match very well in phase and am- 


Fig.4. Phased radial velocities for HD 101930. The planet has 
a minimum mass of 0.30 Mjup and an orbital period of 70.5 
days. 

plitude. Fig|3shows the phased radial velocities for HD 102117 
with the corresponding residuals around the ht, and Table |2l 
gives the best-fit parameters for HD 102117 b. 

5. Discussion and conciusion 

Table |2] summarizes all orbital and physical properties of the 
planets presented in this paper. These new discoveries have 
been made possible thanks to the high RV precision reached 
with the HARPS spectrograph. This precision is demonstrated 
by the very low residuals around the orbital hts, allowing an 
accurate determination of the orbital parameters even with few 
data points, provided the phase coverage is good enough. 

All three stars discussed in this paper are more metal-rich 
than the Sun and therefore conhrm the observational trend 
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Table 5. Radial-velocity values and error bars for HD 102117. 
All data are relative to the solar system barycenter. 


JD-2400000 

RV 

[km s '] 

Uncertainty 
[km s'] 

53017.85639 

49.5717 

0.0014 

53054.80495 

49.5827 

0.0014 

53145.59708 

49.5754 

0.0010 

53147.57142 

49.5843 

0.0009 

53150.67935 

49.5911 

0.0009 

53153.55336 

49.5937 

0.0009 

53156.57259 

49.5871 

0.0009 

53201.49891 

49.5773 

0.0008 

53204.47977 

49.5740 

0.0009 

53207.46689 

49.5783 

0.0013 

53217.52802 

49.5903 

0.0009 

53315.86091 

49.5912 

0.0008 

53400.83082 

49.5935 

0.0010 


HD 102117 harps 
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Fig. 5. Phased radial velocities for HD 102117. The planet has 
a minimum mass of 0.14 Mjup and an orbital period of 20.7 
days. 

showing a clear correlation between planet occurrence and 
metallicity. 

As already mentioned in Sect. 1, these new Saturn-mass 
planets are located at intermediate distances from their parent 
star (0.1-0.5 AU). The existence of planets in this area of the 
mass-period diagram might pose a challenge to the standard 
planet formation theories. The difficulty arises from the fast 
runaway accretion and migration processes that would make 
such planets more massive and come closer to their parent 
stars, thereby becoming hot Jupiters (or hot Saturns). A reliable 
mechanism permitting termination of the accretion process at 
intermediate distances from the star remains to be found, al¬ 
though the evaporation of the disk and the metallicity of the 


central star might play an important role. Among the three plan¬ 
ets presented in this paper, HD 102117 b is certainly the most 
challenging one in that respect, with a mass of ~50 and a 
semi-major axis of 0.15 A U. The recently disc overed Neptune- 
mass planet around fi Ara dSantos et al.l2004^ . orbiting at 0.09 
AU, shows similar properties, but in that case the presence of 
other massive planets in the system might have had an influ¬ 
ence on its final mass. On the contrary, no other massive planets 
seem to exist around HD 102117, at least within ~5 AU. From 
the observational point of view, these new discoveries might 
indicate that a population of Neptune- and Saturn-mass planets 
remains to be discovered below 1 AU. The increasing precision 
of the radial-velocity surveys will help answer this question in 
the near future, thereby providing useful new constraints on 
planet formation theories. 
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